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Abstract 

We first highlight the main differences between second order and higher order linear parabolic 
equations. Then we survey existing results for the latter, in particular by analyzing the behavior of 
the convolution kernels. We illustrate the updated state of art and we suggest several open problems. 
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1 Introduction 

The Cauchy problem in M. n (n > 1) for higher order (m > 2) linear parabolic equations 



' Ut + (-l) m D a {a a ^}u = 



|a|<m 
\$\<m 



u (x, 0) = uq(x) 



m 



m 



has recently attracted some interest, due to its somehow surprising and unexpected properties, strikingly 
different when compared with the corresponding second order parabolic equation, that is, when m = 1. 
The purpose of the present paper is to survey existing results about problem ([1]) and to suggest several 
open problems whose solution would contribute towards the formation of a complete theory. 
Even in the simplest situation when ([1]) becomes the polyharmonic heat equation 



u t + (-A) m u = in 
u (x, 0) = uq(x) in 



important differences appear and many questions are still open. As was first observed by Evgrafov- 
Postnikov [21], the kernels of the heat operators in ([2]) depend on the space dimension, contrary to 
the classical second order heat operator; this apparently harmless fact, already claims a lot of work in 
order to obtain fine qualitative properties of the solution to (|2|). When uq G C° n L°° (M n ), problem 
([2]) admits a unique global in time bounded solution explicitly given by 



u(x,t) = at- n / 2m f Mx-y)fm, n (-M^)dy, (x,t) 



(3) 



1 



where a = a, 



■m,n 



> is a suitable normalization constant and 



l-n 



poo 

/ e~ s2m (r]s) 
Jo 




\n-2)/2(Vs)ds 



(4) 



see |12j . Here and below, J u denotes the v-th. Bessel function. So, not only the kernels / mj „ depend 
on n, but also they are not available in a simple form. Due to the presence of Bessels functions in ([3]), 
the solution to ([5]) exhibits oscillations and this fact has two main consequences. First, the positivity 
preserving property fails; it is in general false that positivity of the initial datum uq yields positivity 
of the solution u. Second, in order to prove global existence or finite time blow-up for corresponding 
semilinear equations, comparison principles cannot be used; for this reason, Galaktionov-Pohozaev |23j 
introduced a new method based on majorizing order-preserving operators which, basically, consists in 
taking the convolution of the initial datum uq with the absolute value of the kernel f m ,n- 

The asymptotic behavior of the solution to the second order heat equation can be described with 
some precision also thanks to the so-called Fokker-Plank equation obtained by exploiting the self- 
similar structure of the fundamental solution. But the Fokker-Plank operator corresponding to ([2]) is 
not self-adjoint if m > 2 and this brings several difficulties to the analysis of its spectral properties; 
these difficulties were partially overcome in a fundamental paper by Egorov-Galaktionov-Kondratiev- 
Pohozaev |19j . However, most of the classical methods usually exploited for the second order heat 
equation do not apply. For instance, any reasonable Lyapunov functional becomes very complicated 
due to the presence of higher order derivatives, too many terms appear and the study of their signs is 
out of reach. Also standard entropy methods fail, due to the change of sign of the kernels f m ,n- the 
second order entropy is J ulogu and cannot be considered because the solution u to ([2]) changes sign 
also for positive data. The sign change of the kernels also forbids to analyze the behavior of suitable 
scaled ratios such as uj ' f m ,n in order to obtain Ornstein-Uhlenbeck-type equations. 

The fact that the functions f m ,n exhibit oscillations also implies that the semigroup associated to 
(PQ) is not Markovian if m > 2; this yields important complications in extending the L 2 theory to an 
LP theory. In the second order case one uses the Markovian properties of the L 2 semigroup to prove 
that it extends to a contraction semigroup in L p . This then leads to heat kernel estimates, a topic 
extensively studied in the past 25 years. For m > 2 and L°° coefficients the situation is reversed: one 
first obtains heat kernel estimates and then applies them in order to develop the L p theory. The heat 
kernel estimates depend essentially upon the validity of the Sobolev embedding H m (M. n ) C C°(R n ), 
hence an important distinction arises depending on the dimension n. This is in contrast to the second 
order case where the theory does not depend on such an embedding. 

The problem of obtaining sharp heat kernel estimates is itself very interesting. To put it into context, 
one needs to go back to short time asymptotic estimates, first proved by Evgrafov-Postnikov [21] for 
constant coefficient equations and later extended by Tintarev [31] for variable smooth coefficients. 
Progress has been made in the past years in obtaining sharp heat kernel bounds, but several important 
questions remain open. 

Further recent results are available for ([2]). In [TUj the positivity preserving property is studied in 
presence of a source f(x,t). In [13] the solvability of the Cauchy problem ([2]) (with m = 2) in presence 
of an irregular datum uq is studied and the presence of a strongly continuous analytic semigroup is 
proved. Finally, we mention that more general linear problems were considered in [30J whereas the 
stability method for higher order equations was studied in [2H Chapter 12]. 

For the above reasons, many natural questions arise. In this paper, we mainly focus our interest on 
the fundamental solution (heat kernel) of (pQ) and on positivity preserving property (ppp from now on). 
As already mentioned, these problems are by now very well understood in the second order case where 
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heat kernels have been extensively studied in very general frameworks, while ppp holds as consequence 
of the positivity of the Gaussian heat kernel (maximum principle). In the higher order case the situation 
is considerably more complicated and it is precisely our purpose to give an updated state of art as well 
as a number of open problems still to be solved in order to reach a satisfactory theory. 

In Section [2] we study various properties of the heat kernel of the general problem ([1]). We avoid any 
local regularity assumptions on the coefficients, and we start with Davies' results |14t \TE[ [T6| [T7] on 
operators with L°° coefficients, omitting reference to earlier work where local regularity assumptions 
were imposed. We then tackle the LP theory, emphasizing the dimensional dependence. We proceed to 
present the short time asymptotics of Evgrafov-Postnikov and Tintarev and sharp heat kernel estimates, 
including results on non-uniformly elliptic operators. In the last part of Section [2] we restrict our 
attention to constant coefficients case, namely equation ([2]). In this simplified situation, especially if 
m = 2, much more can be said on the behavior of the kernels; in particular, we exhibit fine properties 
of their moments. 

In Section [3] we transform ([2]) into a Fokker-Planck-type equation and we recall an important result 
by Egorov-Galaktionov-Kondratiev-Pohozaev [19] about the spectrum of the corresponding (non self- 
adjoint) operator. In Section |4] we determine the behavior of the moments of the solution to the 
Fokker-Planck equation in the fourth order case m = 2. 

In Section [5] we recall the results which describe the way how the ppp may fail and we discuss the 
possibility of finding a limit decay of the datum uq for which ppp may still hold. 

2 Heat kernel estimates 

In this section we survey some properties of the heat kernel of problem (pQ). We first discuss the case 
where the operator has L°° coefficients, then we extend some results to the "singular case" where the 
coefficients are merely assumed to be in L™ c , finally we specialize to the the simplest case of constant 
coefficients for fourth order equations: the biharmonic heat kernel. 

2.1 Semigroup generation 

Problem ([TJ is to be understood in the L 2 -sense, and for this we need to properly define the elliptic 
operator 

(Hu)(x) = (-l) m £ D a {a a ^(x)D^u} 

\a\<m 
\P\<m 

as a self-adjoint operator in L 2 (M n ). For this we start with real- valued functions a a ^(x) = ap jOC (x), 
\a\, |/3| < m, in L°°(M. n ) and we define the quadratic form 

Q(u)= [ V a a)l3 (x)D a uD l3 udx 

\a\<m 
\$\<m 

on Dom(Q) = H m (W n ). Our main ellipticity assumption is that Garding's inequality 

Q(u) > ciHfp^nj - c 2 \\u\\ 2 L2{Rn) , u G H m (R n ), (5) 

is satisfied for some c\,C2 > 0. It then follows that the form Q is closed; the operator H is defined as 
the self-adjoint operator on L 2 (M n ) associated to the quadratic form Q. It is well-known [1, Theorem 
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7.12] that inequality (0) implies that the principal symbol of H satisfies 



2m 



n 



|o|=m 
|/3|=m 



and that the converse implication is true for uniformly continuous coefficients. 

We first consider the question of existence of a heat kernel together with pointwise estimates. The 
heat kernel K(t,x,y) of H is, by definition, the integral kernel of the semigroup e~ Ht , provided such 
kernel exists. Hence it represents the solution u(x, t) of ([!]) in the sense that 



The results depend on whether the order 2m of H exceeds or not the dimension n. 

Theorem 1. [HI Lemma 19], [201 Theorem 1.1], [4, Proposition 28] 

If 2m > n then the semigroup e~ Ht has a continuous integral kernel K(t,x,y). Moreover there exist 
positive constants C{, i = 1,2,3, such that 



for all t £ M + and x,y G M. n . 

One application of this theorem is the extension of the L 2 -theory to L p (M n ). The fact that the 
semigroup e~ Ht is not Markovian makes this problem quite different from the second order case m = 1. 

Theorem 2. [HI Theorems 20 and 21] 

Assume that 2m > n. The semigroup e~ ~ Hz , Rez > 0, extends from L 2 (R n ) n L p (R n ) to a bounded 
holomorphic semigroup T p (z) on L p (R n ) for all 1 < p < oo. Moreover, for 1 < p < oo the semigroup 
T p (z) is strongly continuous and its generator —H p has spectrum which is independent of p. 

In the case 2m < n critical Sobolev embedding into LP spaces appear and the situation is different. 

Theorem 3. [HI Theorem 10] 

Assume that 2m < n. Let p c = 2n/(n — 2m) be the Sobolev exponent and let q c = 2n/(n + 2m) denote 
its conjugate. 

(i) The semigroup e~ Hz extends to a strongly continuous bounded holomorphic semigroup T p (z) on 
L p (M n ) for all q c < p < p c . Moreover the spectrum of the generator —H p of T p (z) is independent of p. 

(ii) Assume that m is even. For p [q c ,Pc] there exists an operator H of the above type for which the 
operator e~ Ht does not extend from L 2 (M n ) n L p (M n ) to a bounded operator on L p (W a ), for any t > 0. 
Ln particular the semigroup e~ Ht does not have an integral kernel satisfying (EJ). 

We note that when m is odd a result analogous to (ii) is valid for elliptic systems [15] . We also note 
that if the coefficients are sufficiently regular then a Gaussian heat kernel estimate is valid without any 
restriction on the dimension; see [16] and references therein for more details. 




K(t,x,y)u (y)dy , V(i,t)ei n xl + . 
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2.2 Short time asymptotic estimates 

In this subsection we make the additional assumption that the coefficients {a a p(x)} are smooth. We 
consider the problem ([1]) and denote by 

\a\=m 
|/3]=m 

the corresponding principal symbol, which satisfies 

for some c > 0. The following notion of strong convexity was first introduced by Evgrafov-Postnikov 
[21j . For a multi-index 7 with I7I = 2m we denote c 2m = (2m)!/(7i! . . . 7 n !)- We define the functions 
6 7 (x), I7I = 2m, by requiring that 

|7|=2m 

Definition. The symbol A(x, £) is strongly convex if the quadratic form 

T(x,v) = Y b a+p(x)v a vp , v = (v a )eC u , 

\a\=m 
|/3]=m 

is positive semi-definite for all x G M n . It is known [211 Section 1] that strong convexity implies that 
the matrix {A^^ (x, is positive definite for all x € W 1 and £ 6 M n \ {0}. 
We first consider operators with constant coefficients so that K(t,x,y) = K(t,x — y, 0). We set 

2m / 7T \ 

<7 m = (2m-l)(2m) 2m_1 sin (^73^) • ( 7 ) 

Theorem 4. [2TJ Theorem 4.1] 

Assume that H is homogeneous of order 2m with constant coefficients and that the symbol A(^) is 
strongly convex. Let 

There exists a positive function S(x) such that for any x £ M. n , x 7^ 0, we have 

2m 

K(t,x,0) = S(x)t ^cosk^ cot — — - -L- o ; +o(l) 



£ 2m-l 
/" 2m 

x exp < -CJ m ; 

I £2m-l 



4m - 2 / 4m - 2 



as i — >• 0. 
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In order not to become too technical we refer to [5T] for the precise definition of S(x); we note however 
that it is positively homogeneous of degree —n(m — l)/(2m — 1) in x G R n \ {0}. 

To extend Theorem H] to the case of variable smooth coefficients we need some elementary notions of 
Finsler geometry. Very roughly, one can say that a Finsler metric is the assignment of a norm at each 
tangent space of a manifold. In our context, extending flSJ) we define 

p(x,0 = max - g \^ /9 , x G R n , £ G R n . 

This defines a Finsler metric on W 1 in the sense that 

p(x,£) = if and only if £ = and p(ai, Af) = | \\p(x, £), A G R. (9) 

In Finsler geometry the definition is typically complemented by 

the matrix { gij } := \ ''-£- ; jy is positive definite for all x G W 1 and £ G R n \ {0}. (10) 

For our purposes we shall not assume (|10p except in Theorem [8] below. We note however that if ©- 
(|10p are valid then the map £ i— >■ p(x,^) is indeed a norm for all x G M n . The length of an absolutely 
continuous path, 7 = j(t), < t < 1, is then defined as 

Z( 7 )= /" p( 7 (t),7(t))dt, (11) 
•/ 

and the Finsler distance between two points 1,1/6 R n is given by 

d(x,y) = inf{Z( 7 ) : 7 has endpoints x and y}. 

Theorem 5. [31, Theorem 1.1] 

Assume that the operator H is homogeneous of order 2m with smooth coefficients and that the prin- 
cipal symbol A(x,£) is strongly convex. Assume further that the matrices {a>a,j3{ x )}\a\=\p\=m and 
{A^ (a;,£)}i<ij<n ar e both positive definite uniformly in x G M n and £ G S n . Then there exist 
functions Vk(t, x,y), k = 0, 1, . . ., such that the following is true: for any x G W 1 there exists 5 > such 
that for < \x — y\ < 5 the following asymptotic expansion is valid as t —> 0: 



fc-f d(x v)2^l I 
K(t,x,y) ~\t*™-iy k (t,x,y)expl-<T m v } . (12) 

27ie functions Vk(t,x,y) oscillate and are bounded and smooth with respect to t. 
Estimate (|12p is meant in the sense that for each N > 1 and for small enough t > there holds 

N 



\ J '~% / s dfa;, y) 

K(t,x,y) - ^i 2 '"- 1 W,£,y)exp ^ -<j„ 



fc=o 



£ 2m— 1 



AT+1- 



< c^t 2m - 1 exp < — o~ r 



d(x,yY- 



2m 



f2m-l 
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2.3 Sharp heat kernel bounds 

We now return to the general framework of operators with L°° coefficients satisfying Garding's in- 
equality ([5]). We assume that 2m > n so that the heat kernel estimate (|6|) is valid and we present 
certain theorems that provide additional information on the constant C2 in ([5]) . The sharpness of these 
estimates is measured by comparison against the short time asymptotics of Theorem [5j 

Theorem 6. [9, Theorem 4.5] 

Let H be an operator of order 2m > n with real-valued coefficients in L°°(]R n ). Assume that the 
principal coefficients {a a> p(x)}\ a \ = \p\ =m satisfy 

a° aj3 VaVp < a a) p(x)vaVf3 < V a° a/3 v a vp, veC, Iff, 

|a|=m |a|=m |a|=m 

|/3|=m |/3|=m \(5\=m 

for some \i > 1, where {0^} is a coefficient matrix for (— A) m . Then for any e > there exists c € 
such that the heat kernel of H satisfies 

{ 2m \ 
ii \ \\ x ~ v\ 2m ~ 1 
-(p(m,n) - e) — h c e t > , 
t2m-l J 

for all t G M+ and x, y G M. n , where 

p(m,n) = (2m - l)(2mr 2m /( 2m - i y /(2m ~ 1) f sin ( — - — V 2m+1 + CyTU - 1 

L V Am — 2) 

and the constant C depends only on m and n. In particular p(m, n) = a m + 0(p — 1) as [i — > 1 + . 

While Theorem [6] provides useful information when H is close to (— A) m , it is clearly not very effective 
when H is an arbitrary elliptic operator. In such a case, the Finsler distance should play a role. Since 
definition (jlip is meaningless when H has measurable coefficients, an alternative definition is required, 
as was the case for second order operators. Denoting by A(x,£) the principal symbol of H we define 

8 = {<f> G C 1 ^") : A(x,V<f>(x)) < 1 for almost all x G R n }. 

For operators with smooth coefficients the Finsler distance d(x, y) is then also given by 

d(x, y) = sup{0(y) - <p(x) : <f> G £} ; (13) 

see [21 Lemma 1.3]. Hence we use (|13p to define the Finsler distance when H has measurable coefficients. 
We note that a simple approximation argument shows that in the definition of £ we could have required 
that <p G C°°(M. n ). Given M > we also define 

£ M = {<Pe C m {R n ) : A{x,X7cf)(x)) < 1 , \V k (p(x)\ < M , a.e. i£K™, 2 < k < m} 

and the Finsler-type distance 

du(x,y) = sup{4>(y) - (p(x) : <p G £m } • (14) 

So doo(x,y) = d(x,y), but for finite M we have dM(x,y) < d(x,y) in general. 
We finally define the following measure of regularity of the principal coefficients of H, 

q A = max dist^^a^, W™' 1 ' 00 ^)). 

\a\=m 



-rn 



In particular q A = if the principal coefficients are uniformly continuous. 
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Theorem 7. Theorem 1] 

Let 2m > n. Assume that the principal symbol A(x,£) is strongly convex. For any M > and e > 
there exists a constant T t m such that the heat kernel of H satisfies 



2m_ 
-1 



\K(t, x, y)\ < T e , M t-^ exp \-(a m - Cq A - e) V ) ^ + T eM t } , (15) 



for all t € M+ and x,y £ 



The constant C in (|15p depends only on m, n and the constants in Garding's inequality ©. In 
relation to the last theorem we mention the following open problems: 



Problem 1. Is the term CqA necessary in (|15p ? Under what assumptions can it be removed? 

Problem 2. Is it possible to replace dM(x,y) by d(x,y) in (|15p ? Under what assumptions can? 

Problem 3. What is the role of strong convexity in the above theorems? What are the best possible 
results if we do not assume the strong convexity? 

Problem 4. For operators with regular coefficients obtain sharp heat kernel estimates when 2m < n. 

A partial answer to Problems [T] and [2] is provided in the next theorem under additional assumptions 
on the principal coefficients. Of course, the questions remain as to what is the best possible result for 
measurable coefficients. The proof of the theorem is geometric and consists in showing that du/d —> 1 
as M — > +oo, uniformly in x and y. 

Theorem 8. [6, Corollary 3] 

Let H be an elliptic operator of order 2m > n whose principal symbol A(x, £) is strongly convex, is 
(jm+i yyftfo reS p ec t to x and satisfies \V^.A(x,^,)\ < c|£| 2m , < k < m + 1. Assume further that the map 



(x,0 A(x,£) 



i 

2rn 



defines a Finsler metric on M. n in the sense that (I?\)- ilQ\) are satisfied. Then the heat kernel of H 
satisfies the estimate 

f 2m \ 

n d( X V ) 2m— 1 

\K(t,x,y)\ < c e t~2^ exp <^ -{a m - e)^^ + c e t }, (16) 



£2m-l 



for any e > and all t G M + and x, y £ 



We next consider singular operators with unbounded coefficients. Let a a ^{x) = ap %a (x), \a\ = \/3\ = 
m, be real-valued functions in L^ c (M n ). We fix s > and assume that the weight a(x) = 1 + |x| 
controls the size of the matrix {a Qj( g} in the sense that 



' 1 a(x)\v\ 2 < ^2 a a,i3(x)v a V/3 < ca{x)\v\ 2 , v G C u , X £ 

\a\=m 

\P\=m 



We consider the elliptic operator 



Hu={-l) m ^ D a {a a , p DPu} 

\a\=m 
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on L 2 (R n ), defined by means of a quadratic form similarly to the uniformly elliptic case; see [5] for 
details. For M > we then define the set 

E M = {<!>£ C m (R n ) : A(x,V(f>(x)) < 1, |V fc 0l < Ma(x)- k/2m , a.e. 2 < k < m} 



and the Finsler-type distance (fT3]l . The weight a(x) induces the weighted L°°-norm ||tt||£,c»rR»») = 
sup R n(|n|/a) and more generally the weighted Sobolev spaces 

W*'°°(R n ) = {u£ Wf£°(R n ) : \V j u(x)\ < ca(x)( 2m -^/ 2m , a.e. x £ R n \ < j < m - 1} 

We set 

g A = max dist L oc (R „ ) (a a , /3 , W™- 1 ' 00 (R n )) . 

\a\=m 
\P\=m 

Theorem 9. Section 2] and Theorem 2.2] 

Assume that n is odd, that < s < 2m — n and £/iai the principal symbol of H is strongly convex. 
Then for any M > and e > there exists a constant T t ^M such that the heat kernel of H satisfies 

{ 2m \ 
-(a m - cq A - 6) dM(Xl ? )2m " + T eM t 1 , (17) 
£2m-l J 

/or a// (€K + and x, y G R n . 

Problem 5. Find out what happens when n is even. 

We end this section presenting a theorem of Dungey [18] for powers of operators. Let (X, d) be a 
metric space and \x be a positive Borel measure on X. Assume that X is of uniform polynomial growth, 
that is there exists c > and D, D* £ N such that the volume V(x, r) of any ball B(x, r) satisfies 

c~ 1 r D < V(x, r) < cr D , if r < 1, 
c~ 1 r D * < V(x,r) < cr D , if r > 1. 

Accordingly let 

Theorem 10. PJ2 Theorem 1] 

iei H be a non-negative self-adjoint operator on L 2 (X,dfi). Assume that the semigroup e~ Ht has an 
integral kernel K(t, x, y) which is continuous in (x, y) for all t £ R+ and satisfies the Gaussian estimate 



\K(t,x,y)\<c e V(ty 



/or any e > and all t £ and x,y £ X. Then for any integer m > 2 the semigroup generated by 
—H m has an integral kernel K m (t,x,y) which satisfies the Gaussian estimate 



\K m {t,x,y)\ < c t V{t)-^ex V \-(o- m -e) d ^' U] 



(2m-l 



for any e > and all t £ R + and x, y £ X . 
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2.4 More on the heat kernel of the biharmonic operator 

In particular situations, much more can be said about the kernels relative to (pQ). In this subsection we 
collect a number of properties related to the heat kernel of the polyharmonic operator (— A) m . All the 
information about the heat kernel of (— A) m is contained in the functions f m ,n, since (cf. ([3])) 

K (t, X,y) = Ot m ^ n t I fm,n ^ £i/2m ) ' 

We specialize to the case m = 2 and we give some hints on how to obtain the corresponding results 
in the higher order case m > 3. For simplicity, we denote f n = /2 n . 
When m = 2, ([2]) becomes the Cauchy problem 

u t + A 2 u = in R n x R + , 

u(x,Q) = u (x) in R n , ^ ' 

whereas the kernels defined in read 

fn(v) = r, l ~ n / e' s M n / 2 J (ri _ 2)/2 M als . (19) 
J o 

These kernels obey the following recurrence formula, see [22J: 

fn(v) = -vfn+2(v) foralln>l. (20) 

Moreover, thanks to Evgrafov-Postnikov [21] (see also [291 (1.10)]), we know that the kernels have 
exponential decay at infinity. More precisely, define the constants 

a = —T7T > K n 



16 ' n (2vr)«/2 . 2 (n-3)/3 ' 

then, in any space dimension n > 1, we have 



= {cos (vW 3 ~)+ 0( V - A/3 )} e"^ 4/3 as r, -+ oo 



(21) 



In [3] one can find the definition of the Gamma function and the power series expansion of the Bessel 
function: 



roo 00 (_-\\k( /<2\2k+i/ 



as well as further properties of T and J u . This allows to obtain the representation of f n through power 
series: 



Theorem 11. [22J Theorem 2.1] 

For any integer j > 1, we have 



r[i±i 



Mi) = E tm«m is ' <22) 
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For any nonnegative integer j, we have 



h (fc + j)!r(2*±fi±i 

-l)k V * 



Wo) = £ E <-d' ' ; ■ < 23 > 

In particular, f n (0) > for all n and 

i ^ . r fasti) ,„ ... i » . , r(i±i- 



AW = ^E • aw = jE (-')' |^ 

v A;=0 V ; fc=0 1 J 



v 2k 



Using the properties of the Bessel functions, the following third order ODE for the function f n was 
derived in [22\ Theorem 2.2] for any integer n > 1: 



or, equivalently, 



fniv) + ^fn(v) ~ ^fn(.V) ~ |/nW = (24) 

(Af n )'( V ) = ^f n ( V ) . (25) 

According to (|2ip the kernel f n (jj)i an d hence the biharmonic heat kernel, has infinitely many sign 
changes as rj — > oo, see also previous work by Bernstein [11] when n = 1. We refer to [2H] for further 
(minor) properties concerning the behavior of the kernels at some special points. 

We now rescale the kernel f n and define the function: 

/•oo 

v 0O (y) = 2 n / 2 a n f n (V2\y\) = 2 ( - n + 2 y*a n \y\ 1 - n / 2 e~ s * s n ' 2 J [n _ 2)/2 {V2\y\ s) ds \/y G W 1 (26) 

Jo 

where a n is given by 



a, 



l =u n r n - 1 f n {r)dr= / f n (\x\) dx 



o 



here uo n denotes the surface measure of the ra-dimensional unit ball (so that oj\ = 2). Note that 
j^nVoo(y)dy = 1. Although the functions Voo and f n are strictly related we maintain the double 
notation since, in our setting, they play quite different roles; the former is a stationary solution to (|34p 
below, the latter is the biharmonic heat kernel. We aim to study the moments of the function Voq 
defined in (|26p . The prototype monomial in R n is given by 



Pi(v) = y e -= n y l i for 1 = ~>40 G N " ( 2? ) 

i=l 

and its degree is \i\ = Yli^i- Then we define the P^-moment of Voo by 

M P( := ( Pe(y)v 00 (y)dy (28) 

and we have 

Theorem 12. [251 Theorem 2] 

For any £ = (£i, £ n ) E N n the following facts hold: 
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1. M A 2 Pe = -\£\M Pe , 

2. if \(.\ 4N or if at least one of the 1% 's is odd, then Mp t = 0, 

3. if \l\ G 8N and all the li 's are even, then ■M.p l > 0, 

4- if \l\ S 8N + 4 and all the li 's are even, then Mp e < 0. 

We have so far considered moments having polynomials of y as weights; we now consider powers of 
\y\ which are polynomials only for even integer powers. For any b > —n we define 

M b := [ {yfvMdy . (29) 

Note that for b > — n the above integral is finite since \y\ b v oc (y) ~ foo(0) \y\ b as y — > and has 
exponential decay at infinity according to (|2ip and (I26p . If Pe(y) = \yf for some I G 2N, then 
coincides with Aip t as defined in (|28p . We are again interested in the sign of these moments. The 
following result holds: 

Theorem 13. [251 Theorem 4] 

Assume that n > 1 and that b > —n. Then 

oo 

M h > , for all b £ (-n, 2) [J ( |J (8fc + 6, 8A; + 10)) , 

fc=0 

7W 6 = 0, for all b G 4N + 2 , 

oo 

Xfe < , /or a// |J (8/c + 2, 8A: + 6) . 

When b £ (— n,0], Theorem 1131 was first proved in [22, Proposition 3.2]. Theorems 1121 and 1131 give 
further information about the sign-changing properties of the kernels f n (recall (|26p ). and they better 
describe how these infinitely many sign changes occur. They also show that the sign of the moments 
of f n do not depend on n. 

We conclude this section by explaining how the just described properties of the biharmonic heat 
kernels can possibly be extended to higher order polyharmonic kernels. First of all, we recall that |21| 
Theorem 4.1] (see also |29^ (1.10)]) gives the following generalization to (I2ip in any space dimension 
n > 1: 

fm,n(v) = ^m-m^-l) H {^'^ ~ W) + O )} e"^^ (30) 

as rj — y oo for some (explicit) positive constants K m>n and b m ^ n depending on m and n, and some 
(explicit) positive constants a m and a m depending only on m. 
Next, we suggest the following 

Problem 6. Determine a power series representation of the kind of Theorem [TT] for the higher order 
kernels / m>n (m > 3) defined in @. 
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To this end, by arguing as in [22J and using Section 4.62], it may be useful to notice that (|2Uj) still 
holds, independently of to. Moreover, the following (2m — l)-order differential equation holds: 



(A^/m.n)'^) = t^vUAv) for all n > 1. (31) 

It is straightforward that (|3ip coincides with (|25p if to = 2, whereas it reduces to /'(??) = —^vfiv) 
whenever to = 1 (recall that in the latter case, the kernel / is independent of n). 
With these two identities, one obtains results similar to Theorem 1131 In particular, one has 

Cm,n,/3 '■= w n / r^™ -1 "/ 3 / mn (r/) dr] > for all integers n > 1 and all /3 £ [0,n) (32) 

where w n denotes the measure of the unit ball in M. n . The proof of (|32p can be obtained following the 
same lines as [221 Proposition 3.2], see [25] . 

Problem 7. Prove the full extension of Theorem 1131 to the case of general to > 2. What are the signs 
of C m ^ n fi for all /3 € (— oo,n)? How do they depend on to? 



3 The Fokker-Planck equation 

In some situations it is convenient to transform ([2]) into a Fokker-Planck-type equation. Let 

R(t) := (2mt + l) 1/2m 

so that R(t) 2m - 1 R'(t) = 1. Also put 



u(x, t) := R(t)~ n v [j^^g R{t)\ . (33) 

Then take r = logi?(t) and y = x/R(t). Some lengthy but straightforward computations show that 
v = v(y, t) solves 

' Vt + C v = inl"xl + , 
v(y,0)=u (y) inW 1 , [6 ' 

where 

Cv := {-A) m v- V- (yv) . (35) 



We recall here some properties of the operator C defined in (J35J). The most relevant one is that, 
contrary to the second order heat equation, the operator C is not self-adjoint: we refer to [191 Section 
3] for some properties of the adjoint operator C * . Let o m > be as in ([7|)- ([30|) and, for any a £ [0, <r m ), 
consider the function 

p a (x) = e a ^ ml{2m - X \ xeR n (36) 

so that, in particular, p a = 1 if a = 0. For any such function p a consider the space L 2 (M n ), the weighted 
L 2 -space endowed with the scalar product and norm 

(u,v) L 2(^ n) = p a (x)u(x)v(x)dx , |Ml£2(R«) = («. «)l^(H») ■ ( 37 ) 
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Clearly, if a = we have L^(M n ) = L 2 (R n ). Together with the space L^(]R n ), we consider the weighted 
Sobolev space H 2m (W l ) endowed with the scalar product 

(u,v) H 2 m{R n) = / p a [x) D a u(x) D a v(x)dx . 

|a|<2m 

By |19} Proposition 2.1] we know that £ is a bounded linear operator from .£f 2m (IR n ) onto L 2 (K n ). 

We now wish to characterize the spectrum of £ . In particular, the kernel of £ is nontrivial; any 
function in the kernel is a stationary solution to (134p , As for (|26p . we rescale the kernels f m ,n by setting 

foo(y) = C m ,n/ m ,„((2m) 1/2m |y|), Vy G R n (38) 

where C m ^ n > is a normalization constant chosen in such a way that j^ n v 00 (y) dy = 1; note that 
Vqo G 5, where 5 is the space of smooth fast decaying functions: 

S := {w G C°°(K n ) : |x| a £) a w(x) -> as |x| -> oo for all a > , a G N n } . (39) 

In fact, there exists a unique stationary solution to (j34"j) which belongs to S: 

Theorem 14. [191 Theorem 2.1] 

Up to a multiplication by a constant, there exists a unique nontrivial stationary solution to |g^[ ) which 
belongs to S. This solution v is radially symmetric and, if we further assume that J* Kn v(y) dy = 1, it 
is explicitly given by Voo in [38\). 

Moreover, the spectrum of £ coincides with the set of nonnegative integers, o~(£) = N. Each eigen- 
value A G a(C ) has finite multiplicity and the corresponding eigenf unctions are given by 

D a v OQ for \a\ = A G N . 

The set of eigenf unctions is complete in L^(M n ) for any a G [0, a m ). 

This fundamental (and elegant) result certainly deserves more investigation. Consider the (normal- 
ized) projection operator P a defined by 

P a w:= ( [ p a w Voo dx) ^ for all w G L 2 a (R n ) . (40) 

We recall two problems suggested in [25] . 

Problem 8. Prove the generalized Poincare-type inequality 

[|u - Paullia^n) < (u, C u) L 2 {Rn) for all u G ^ m (IR n ) . 

Although from Theorem Q3] we know that the least nontrivial eigenvalue of £ is 1, since C is not 
self-adjoint the above inequality is by far nontrivial. In particular, prove (or disprove) the following: 

(u,Cu) L 2(wn-) = / p a (x) u(x) £ u(x) dx > Hull^cjjn) for all u G [ker£] _L . 

Problem 9. Determine the convergence rate in TP (for 1 < p < oo) of the solution to (|34p towards its 
projection onto the kernel, that is, onto the space spanned by Voq. 
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4 Asymptotic behavior of the solution 

In this section we shed some light on the long-time behavior of solutions to (|18p . The asymptotic 
behavior is better seen in the Fokker-Planck equation. When m = 2, (|34p reads 

v T + A 2 v -V ■ (yv)v = inK"xK + , 

v(y,0) = u (y) in M" . ^ ' 

We now study the moments of the solution v to ([4"T]) . Let S be as in ([3"9"]) , let uq £ S and consider the 
solution v to ([31]). Let be as in ([27]) and consider the (time-dependent) map 

Mp e , Uo (r):= P e (y)v(y,r)dy = y e v(y,r)dy. 

Let Vqo be as in ([26]) and let M.p e be as in ([28]) . We have 
Theorem 15. [25] Theorem 3] 

Assume that uq £ S is normalized in such a way that 



uo{y)dy= v oa (y)dy = l (42) 

and let v denote the solution to i| ), For any r > 0, i/ie following facts hold: 

(i) M^ >uo (r) = - M A2p ^ tt0 (r) - |i| M Pf , Uo (r) /or a// £ G N» 
(m) Mp, w (r) = e^ |r j Rn Pi(x) u (x) dx for all \l\ < 3, 
{Hi) lim^oo M PltU0 (T) =M Pe for all I G W. 

By combining Theorems 1121 and 1151 we infer 

Corollary 1. Assume that uq G S is normalized in such a way that \J^2^ holds and let v denote the 
solution to i41\ ). Then 

{= if \£\ G" 4N or if at least one of the £{ 's is odd, 
> if \£\ G 8N and all the li 's are even, 
< if \£\ G 8N + 4 and all the li ! s are even. 

In the particular case where \l\ = 2k and Pe(y) = \y\ 2k we may give a simple characterization of the 
moments of a solution to (|4ip . Consider a solution v to (|4ip with initial data uq G S. For all 6 > let 
A^b be as in ([29]) and put 

M b , Uo (T) := I \y\ b v{y,r)dy . 

We then have 

Theorem 16. [25] Theorem 5] 

Assume that uq G S is normalized in such a way that fl^jgp holds and let v denote the solution to |^ j[ ). 
Then for any k G N, k > 2, the above defined functions satisfy the following ODE 

M2k,u (r) + 2k M 2k>Uo (r) = -2k (2k - 2) (2k + n - 2) (2k + n - 4) M 2k -4,u (r) . (43) 



15 



Moreover, for any k £ N, we have 

lim M 2fciU0 (r) =M 2k (44) 



and the following explicit representation 

k 



M 2k ,u (T) = Y; a j e ~ 2jT > ( 45 ) 
j=0 



where a$ = M. 2 k an d 



k-2 k-2 

a j 



(i) a\ = M 2fc , uo (0) + 2k (k - 1) (2k + n - 2) (2k + n - 4) V . 

i=o J 

(ii) = ifk>l, 

k 2fc (fc - 1) (2/c + n - 2) (2fc + n - 4) fc _ 2 

(m) aj = : a - tj k > 2 and j = 0, fe — 2 

k — J 

In (i) we use the convention that ^ji=o = ^ ^ — 
Formula (1451) shows, for instance, that 



M),« ( r )= / u o(y) d V > Af 2 , Uo (r) = e 2t / |y| 2 n (y)dy, 

JR™ JR™ 

M^ U0 ( T ) = -2n(n + 2) j u (y) dy + e~ Ar [ \\y\ 4 + 2n (n + 2)]u (y) dy , 

Jk» jr™ l j 

M 6 , uo (r) = -6(n + 4)(n + 2)e" 2T / |y| 2 u (y)d2/ 

+e" 6T ( / |y| 6 «o(y) dy + 6 (n + 4) (n + 2) / |y| 2 u (y) dy 



If 2N (so that is not a polynomial) we may still define the map Mb jUQ and, for all b E [4, 00), 
we obtain 

K,u (t) + 6 M Mo (r) = - b (b - 2) (b + n - 2) (b + n - 4) M 6 _ 4)Uo (r) . 

Note that Theorems 1151 and 1161 also hold in a weaker form if uq G L x (IR n ) and |?/| a uo € L 1 (M n ) for 
some a > 4. In this case, the statements hold true under the additional restriction that \£\ < a. In 
particular, we have the following 

Corollary 2. Assume that (l + |y| 4 )uo € L 1 (M") and that j^[ ) holds. If v denotes the solution to J^-?[ ), 



lim / |y| 4 t>(y,r)dy = 7W 4 < . 
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5 Positivity preserving property 



Contrary to the second order heat equation, no general positivity preserving property (ppp in the 
sequel) holds for the Cauchy problem ([!]). By ppp, we mean here that positivity of the initial datum 
no implies positivity (in space and time) for the solution u = u(x,t) of ([I]); this is of course equivalent 
to the kernel K(t, x, y) being non-negative. 

Nevertheless, by exploiting the properties of the kernels, some restricted and somehow hidden versions 
of ppp can be observed for the fourth order parabolic equation 

u t + A 2 n = in R"xl + , 

u(x,0) = uq(x) in R n , ^ ' 

where n > 1 and uq G C° n L°° (W 1 ). In this section we recall several weakened versions of ppp for 
the problem (|46p . We start however with a theorem about the general problem (pQ) which provides 
quantative information on the positivity of the heat kernel near the diagonal {x = y}. 

Theorem 17. [17^ Theorem 6] 

Let H be an homogeneous elliptic operator of order 2m > n acting on L 2 (W l ). There exists constants 
ci, C2 > such that 

K(t,x,y) > cit~^ (47) 
for all t 6 1R + and x, y 6 M. n such that \x — y\ 2m < c^t- 

Theorem [T71 states that the solution u = u(x,t) to (P) when Uo(x) = 5{ x=z y (the Dirac delta distribu- 
tion at some z E W 1 ) satisfies u(x,t) > whenever \x — z\ 2m < c^t. Therefore, one expects that if the 
mass of no is "concentrated" in some small region of M. n then ppp holds, at least in some part of that 
region. This can be made precise for the simplified problem fj46f) on which we focus our attention for 
the rest of this section. 

Theorem 18. [271 Theorem 1] 

Assume that ^ uq > is continuous and has compact support in W 1 . Let u = u(x,t) denote the 
corresponding bounded strong solution of j^g[ ). Then, 

(i) for any compact set K C W 1 there exists Tk = Tk{uq) > such that u(x, t) > for all x G K and 
t > T K ; 

(ii) there exists r = t(uq) > such that for all t > r there exists xt £ W 1 such that u(xt,t) < 0. 

The trivial example no = 1 shows that, at least for statement (ii), the compact support assumption 
cannot be dropped. By Theorem 1 181 we see that negativity for exists in general and goes to infinity. 
Fine results concerning the validity of the eventual positivity property in presence of a source, may be 
found in [ID]. 

It appears instructive to combine Theorem 1181 with the following energy conservation laws obtained 
in \25\ Corollary 1]: let no 6 L 1 (M n ) and let u be the solution to (jUJ); then, for all t > we have 



u(x,t)dx= I uo(x)dx, (48) 
f u(x, t) 2 dx = - 2 f \Au(x,t)\ 2 dx . (49) 
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Denote by u + = max{it, 0} and u~ = — min{n, 0} the positive and negative parts of a function u, so 
that u = u + — u~ . Theorem 1181 states that if no £ C°(IR n ) has compact support and ^ uq > in R n , 
then u~(x,t) ^ for all t > 0. Moreover, (|48|) states that the map 

1 1-> / u(x,t)dx (t > 0) 

JR™ 

is constant and equals a strictly positive number. Hence, 

/ u~(x,t)dx> / u~(x,0)dx = for all t > , 

JR n il" 



n + (x, i) (fx > / n + (x,0)ci3; 



/ uq(x) dx for all t > : 



here we use redundant notations (u + (x,0) = Uq(x) = uq(x) and u~(x,0) = Uq{x) = 0) in order to 
emphasize the strict inequalities between the mass of the positive (respectively, negative) part of the 
solution u = u(x,t) and the the mass of the positive (respectively, negative) part of initial datum no- 
On the other hand, (H2J) states that 

t ^ / u(x,t) 2 dx (r > 0) 

decreases and, in particular, that 

u + (x,t) 2 dx< / uq(x) 2 dx = I u^(x) 2 dx (t > 0) . 



Summarizing, the L 2 -norm of the positive part of the solution u is smaller than the L 2 -norm of the 
positive part of the initial datum no, whereas the L 1 -norm of the positive part of the solution u is 
larger than the L 1 -norm of the positive part of the initial datum uq. 

Problem 10. Prove the counterpart of Theorem [181 for ([2]) (for any m > 2) when ^ no > is 
continuous and has compact support in W 1 . 

Next, we consider initial data no which are not compactly supported and which display a given decay 
behavior as \x\ — > oo. We fix some arbitrary /3 > and consider the functional set 

Ca := {g £ C°(IR n ;M + ) : 5(0) > 0, g(x) = o(\xf) as \x\ -> 00} . 

In a suitable class of initial data, a positivity result for the linear Cauchy problem (|46p holds: 

Theorem 19. [221 Theorem 1.1] 
Let (5 > and let g G Cp. Let 

g(x) + \x\P 

Let u = u(x,t) be the corresponding solution of [(6\ ) an d K Cl" be a compact set. 
(i) If j3 < n, then there exists C n> p > such that 

lim ^ /4 n(x, t) = C n s, 

t— >+oo 
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uniformly with respect to x G K . 
(ii) If /3 >n and g(x) = 1, then there exists D n ^ > such that 

lim t n /\\ogt)- 1 u{x,t) = D n>n iff3 = n 
*^+°° ~ (51) 

lim t n/i u(x, t) = D n[S if/3>n, K ' 

t— >-\-oo 

uniformly with respect to x G K . 

Problem 11. By using (|32p . prove the counterpart of Theorem 1191 for ([2]) when no is as in (|50|) . 

The constants C n> p and D n fi m Theorem [19] do not depend on K. What does depend on K is the 
"speed of convergence", namely how fast t^^u(x,t) — C n $ converges to (and similarly for D n ^). 
Let us also mention that if /3 > n, then for any g G Cp (not necessarily constant) one still has that 
lim^+oo t^^u(x, t) = +oo uniformly with respect to x G K. 

Remark 1. The quantitative positivity result of Theorem 1191 provides strong enough information to be 
applied also to semilinear problems, see [22j[26]. At a first glance, this appears somehow unexpected, 
since the techniques connected with the proof of Theorem 1191 seem to be purely linear. 

Theorem [19] does not clarify whether the eventual positivity for solutions of ([46 p is global or only 
local. Theorem [18] suggests that negativity for the solution of (I46D always exists and shifts to infinity, 
provided j3 is sufficiently large. 

Problem 12. Prove Theorem [181 (ii) for any no as in (|50p for (3 large enough. 

On the other hand, if no = 1 then the solution of (|46p is n (x, t) = 1. This trivial example shows that 
if P = 0, presumably one has global eventual positivity for ([46]). At least in the case n = 1, this is also 
true if f3 is positive but sufficiently small: 

Theorem 20. [22] Proposition A. 6] 

We assume that n = 1 and uo(x) = \x\~@ . For (3 > sufficiently small, the corresponding solution of 
(4^\ ) given by 

u(x,t) = a n [ ^ n J,} | fl dz 

is positive in R x M + . 

Problem 13. Prove Theorem 1201 in any space dimension n > 1. 
By combining (|48p with Corollary [2] and with Theorem [20] we obtain 

Corollary 3. Assume that uq > a.e. in W 1 . 

(i) If (1 + |x| 4 )no G L 1 (M n ), then the solution u to M8\) changes sign. 

(ii) If n = 1, there exists (3q > such that if (3 G (0,/3o) and uq(x) = \x\~@ , then the solution u to 
U8\) is a.e. positive inM x R + . 

Corollary [3] can be interpreted as follows. From Theorem 1181 we know that solutions n to (118p with 
compactly supported nonnegative initial data no display the eventual local positivity property, that 
is, u(x,t) becomes eventually positive on any compact subset of M. n but it is always strictly negative 
somewhere in a neighborhood of |x| = oo. This happens because the biharmonic heat kernels exhibit 
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oscillations and, outside the support of uq, they "push below zero" the initial datum. The same happens 
if uq > but uq is "very close to zero", see statement (i). On the other hand, if uq > and no is 
"far away from zero" then the kernels do not have enough negative strength to push the solution below 
zero, see statement (ii). The trivial case uq = 1 (which is a stationary solution to (fl8j) !) well explains 
this situation. 

Finally, the following result shows that in general, we cannot expect neither global positivity nor 
uniform bounds for eventual positivity. 

Theorem 21. [221 Theorem 1.2] 

Let (3 £ (0, n) . For any T > 1 there exists g £ Cp such that if 

«o(z) = , , , g 
g{x) + \x\ p 

then, the corresponding solution u = u(x, t) of |^6p satisfies u (xt, T) < /or some xt £ R n . 
Problem 14. Extend Theorems |2"01 and |2"T1 to (J2J) for any to > 2. 
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